Actin based propulsion: Intriguing interplay between material properties and growth processes 
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Eukaryotic cells and intracellular pathogens such as bacteria or viruses utilize the actin polymerization ma- 
chinery to propel themselves forward. Thereby, the onset of motion and choice of direction may be the result of 
a spontaneous symmetry-breaking or might be triggered by external signals and preexisting asymmetries, e.g. 
through a previous septation in bacteria. 

Although very complex, a key feature of cellular motility is the ability of actin to form dense polymeric 
networks, whose microstructure is tightly regulated by the cell. These polar actin networks produce the forces 
necessary for propulsion but may also be at the origin of a spontaneous symmetry-breaking. 

Understanding the exact role of actin dynamics in cell motility requires multiscale approaches which cap- 
ture at the same time the polymer network structure and dynamics on the scale of a few nanometers and the 
macroscopic distribution of elastic stresses on the scale of the whole cell. 

In this chapter we review a selection of theories on how mechanical material properties and growth processes 
interact to induce the onset of actin based motion. 

PACS numbers: 87.17.Jj 87.15.Rn 82.39.-k 82.35.Pq 62.40.+i 81.40.Jj 



I. INTRODUCTION 

Most living cells are able to perform a directed motion, ei- 
ther by swimming in a liquid environment, by crawling on 
a solid support or by squeezing through a three-dimensional 
matrix of fibers. 

The speed of swimming bacteria can reach up to 
100 |jm s~^, whereas eukaryotic cell crawling can be as fast as 
1 fims^^ (Ref. |fl?|l and references therein). Given their size 
and speed, the motion of single cells is governed by viscous 
forces, not inertia, i.e. the Reynolds number Re <C 1. 

Unicellular organisms move in search of a food or light 
source. Intracellular pathogens like bacteria or viruses spread 
by exiting their host cell and entering a neighboring cell. 
Other simple organisms like the slime mold Dictyostelium 
discoideum migrate under unfavorable conditions, e.g. star- 
vation, towards an aggregation center to form a multicellular 
organism. 

Typically, in a multicellular organism not all cells are motile 
all the time but they can be mobilized by the appropriate stim- 
uli. For example, the ability to move plays a crucial role dur- 
ing embryonic development, in wound healing and in the im- 
mune response. Additionally, cellular motility is a prerequi- 
site for metastasis formation during cancer development 



The biological realizations to produce a propulsive force 
are diverse. Most swimming cells, e.g. sperm cells or the bac- 
terium Escherichia coll, use one or multiple beating flagella', 
respectively. In contrast, crawling cells and some intracellular 
pathogens advance by actin polymerization. 

In this article we will present recent experiments and 
concepts to understand the latter mechanism, the production 
of forces in the advancing edge of crawling cells or for 
the propulsion of intracellular organelles, which is also of 
fundamental interest for the medical and engineering sciences. 

The foundation of the research of cell motility as a dis- 
tinct discipline were laid in the 1970's by the group around 
Michael Abercrombie iQl]. He was the first to divide the mo- 
tion of fibroblasts^ into three phases. Extension - Adhesion - 
Contraction, which form the dogma of cellular motion as it 
is recognized today. In this mechanism, a slow movement 
(~l|jmmin^^) is generated by the extension of flat mem- 
brane sheets, lamellipodia, into the direction of movement. 
The advancement of the membrane is accompanied by the for- 
mation of focal adhesions, contacts between the substratum, 
the cell membrane and actin stress fibers. Finally the cell rear 
retracts, accompanied by a deadhesion of the membrane from 
the substratum. 
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FIG. 1: Actin polymerizes into an elastic filament network (shaded 
in grey) at the leading edge of motile eukaryotic cells or at the outer 
surface of pathogens and organelles and induces cellular motion. The 
direction of motion is indicated by the black arrows. 



The molecular basis and minimal ingredients of actin poly- 
merization forces are now very well understood thanks to 
biomimetic experiments in the group of Marie-France Carlier 
ll33ll . However, the theoretical understanding of the physical 
mechanisms of polymerization process are still a matter of de- 
bate and shall be the main topic of this review. 

In the first introductory part we will outline the biochemical 
basis of actin polymerization and present a selection of exper- 
imentally observed phenomena. Here we will shortly present 
the processes taking place at the leading edge of locomoting 
cells and then we will mainly rely on biomimetic experiments. 
In a second part we will present some theoretical concepts to 
interpret and understand the existing experiments. The list of 
presented models is by no means complete and rather repre- 
sents a selection of the leading ideas. Here we will focus on 
the class of Brownian ratchet models and macroscopic models 
on symmetry-breaking. 



However, there exist variations of this dogma. Fish kera- 
tocytes-' perform a rapid continuous motion (~10nmmin~^) 
with a constant shape. They almost seem to glide over the 
surface and form only transient focal contacts with the sub- 
stratum with a much shorter lifetime than the focal adhesions 
formed in fibroblasts 

Another variant is the rapid motion (~10|jmmin^^) of the 
slime mold Dictyostelium discoideum, which moves in an 
amoeboid fashion. During this amoeboid motion only non- 
specific contacts with the substratum are formed and actin 
stress fibers are absent ifTsll . 

Despite the various mechanisms, cell motion requires first 
of all the self-organization of the cell into an advancing and re- 
ceding edge. This manifests itself by different molecular con- 
centrations or activation levels of enzymes at the two poles. 
The polarization can be guided by external signals, e.g. chem- 
ical gradients or a variation in the mechanical properties of 
the support, but it might also arise in a homogeneous environ- 
ment, after a transient mechanical perturbation of a stationary 
symmetric cell ||60ll . 

The three processes, extension, adhesion and contraction 
are coupled by sophisticated and complex mechanisms. How- 
ever, it seems as if the process of front extension relies on a 
completely different machinery than that of adhesion and con- 
traction and can be studied separately. 

It had been known for a long time that on a cellular level, 
forces can be generated on the basis of muscle like proteins, 
i.e. actin and myosin, which is indeed responsible for the con- 
traction of the cell rear IztIi . However, more recently it was 
discovered that the protrusions at the leading edge as well as 
the motion of intracellular organelles f2% or pathogens, like 
the bacterium Listeria monocytogenes 15311 or the Vaccinia 
virus |J3il , can be associated with the so-called actin poly- 
merization machinery lf52ll as shown schematically in Fig. [T] 



II. EXPERIMENTAL OBSERVATIONS 

A. Biochemistry of actin polymerization and organization of 
the leading edge of advancing cells 

Actin is a sma ll g lobular protein of 42 kDa present in all 
eukaryotic cells IdIII . Under physiological conditions actin 
monomers (G-actin) polymerize into long helical filaments 
(F-actin). In a living organism these polymerization or de- 
polymerization processes are tightly regulated. The literature 
on the dynamics of actin and the proteins which interact with 
actin is vast. As an introduction we would like to refer the 
reader to the comprehensive biochemical reviews by Pollard 
et. al. SlilandRafelskietal. iH or the books by Bray 101 
and Howard ll26ll . Due to the limited scope of this bookchapter 
we will just present the basic phenomena and common termi- 
nology associated with the actin polymerization machinery, 
which will allow the reader to understand the pertinent ques- 
tions and concepts. 



1. Actin polymerization in vitro 

Under physiological conditions, i.e. at an ionic strength of 
^lOOmM, monomeric actin polymerizes spontaneously into 
filaments. The filament growth process typically starts with 
a nucleation process, since actin dimers and trimers are un- 
stable. Shortening or elongation of existing filaments occurs 
predominately via subunit addition or subtraction at the fil- 
ament ends and not via filament breaking or annealing pro- 
cesses. Actin monomers at a concentration c may bind to a 
filament end with a rate ^ k+c and dissociate with a rate 
~ fc_. In a stationary situation, i.e. zero net growth of the 
filament, one can write 

= fc+c-fc_. (1) 



^ fish scale cells 



The concentration = k_/k^ associated with this equilib- 
rium is called critical concentration. 
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FIG. 2: Actin subunits treadmill through a filament. ATP-actin poly- 
merizes at the barbed end, ATP hydrolyzes and ADP-actin depoly- 
merizes at the pointed end. In the solution ATP is exchanged for 
ADP at the monomers. 

G-actin has a structural polarity jzsll . The same polarity 
is also found in actin filaments and can be visualized as an 
arrowhead pattern by decoration of filaments with myosin"*. 
From this arrowhead pattern the two filament ends are referred 
to as pointed and barbed, respectively, as shown in Fig.|2] 

In the presence of Mg-ATP the structural difference trans- 
lates into a difference of the critical concentrations and rate 
constants between the barbed and pointed ends and causes a 
treadmilling of subunits through the filament. 

Briefly, this phenomenon can be explained as follows (for 
details see Ref. |46]). Most actin monomers are bound to ATP 
(typically Mg-ATP). The critical concentration for this ATP 
bound species is about 6 times lower for the barbed end than 
for the pointed end. For ADP actin the critical concentrations 
are about the same for both ends but about ten times higher 
than for ATP actin at the barbed end. Therefore, in the steady 
state the ATP actin concentration is above the critical concen- 
tration of the barbed end and below the critical concentration 
of the pointed end. Polymerized actin subunits are still bound 
to ATP but in the course of time ATP hydrolyzes irreversibly 
into ADP-nPi and, later on, the anorganic phosphate P^ disso- 
ciates from the filament with a half-time of several minutes. 
Consequently, ATP actin polymerizes at the barbed end, trav- 
els along the filament whereby ATP is hydrolyzed and finally 
ADP actin depolymerizes at the pointed end. Fig. |2] shows 
schematically the main key processes of the treadmilling cy- 
cle. 

The process of irreversible ATP hydrolysis and P; release 
(and subsequently the exchange ADP^ ATP at the actin 
monomer in solution) keeps the system out of equilibrium and 
allows for a constant flux of monomers through the filament 
at constant filament length, which is called treadmilling and 
forms the base of cellular motiUty. 



2. Proteins that regulate actin polymerization in vivo 

In the living cell the actin polymerization machinery is 
tightly regulated by signalling processes. Many different fac- 
tors interact with actin or participate in the regulation of this 



Myosin molecules bind to each actin subunit in the filament in an oriented 
fashion, leading to a typical pattern along the filament, which looks like a 
series of arrowheads in an electron micrograph. For details see Ref. 0|. 



polymerization machinery. In the following paragraphs we 
will only discuss the relevant factors which determine the 
actin architecture of the leading edge of advancing cells or 
the actin comets formed by intracellular pathogens and which 
seem to be crucial to generate motion. 

First of all, most actin monomers are bound to so-called 
monomer binding proteins, e.g. Thymosin-/34 and Profilin, 
and are thus not able to nucleate new filaments. But Profilin- 
ATP-actin complexes elongate existing filaments at the barbed 
end nearly as efficiently as ATP-actin. 

Electron micrographs of the leading edge show a dense net- 
work of actin filaments linked to each other by Y junctions, 
where a new filament grows from an existing filament at a 
70° angle. The distance between two crosslinks is of the or- 
der '^20-30nm (in comparison, the persistence length of actin 
filaments is 15 |am ll42ll ). This Y junction is initiated by the in- 
teraction of the Wasp/Scar protein and the Arp2/3 complex, 
whereby Wasp/Scar is a membrane bound protein which in- 
cites the binding of the Arp2/3 complex to an existing fila- 
ment, which in turn then serves as a nucleation point for a 
new filament. 

Free barbed filament ends are quickly covered by so called 
capping proteins, thus limiting filament elongation to a zone 
near the plasma membrane. At the pointed filament end de- 
polymerization takes place. This process can be accelerated 
by a protein complex called ADF/cofilin, which is able to 
sever old filaments containing ADP actin subunits. These fil- 
ament fragments are then depolymerizing rapidly. Besides 
Arp2/3 there exist other types of filament nucleators, called 
formins. In the presence of profilin and ADF/cofilin they seem 
to favor the growth and bundling of several actin filaments into 
cables IJsll . present in spike-like membrane extensions, called 
filopodia. In this chapter we will limit ourselves to actin net- 
works produced by the Arp2/3 complex. 

The above outlined process of polymerization at the mem- 
brane/actin network interface and depolymerization far away 
from the membrane provides the mechanism for pushing the 
membrane into the direction of movement. However, little is 
known about the mechanical properties of such filament net- 
works and recent experiments indicate that the loading history 
determines the growth velocity of the network ll43ll . 

Certain bacteria invade other living cells and hijack the 
actin polymerization machinery of their host cells to propel 
themselves forward. They carry a protein on their outer sur- 
face (e.g. ActA for Listeria monocytogenes), which adopts 
the same function as the Wasp/Scar protein in the membrane 
of eukaryotes: it triggers the polymerization of an Arp2/3 
crosslinked network at the outer bacterial surface. This actin 
network typically develops asymmetrically only at one side of 
the bacterium thus pushing the bacterium in the other direc- 
tion. 

A similar mechanism might also be responsible for the mo- 
tion of endocytic vesicles in living cells (see [29] and refer- 
ences therein). 

To summarize the growth processes and the resulting archi- 
tecture of the actin system at the leading edge: close to the 
cell membrane containing an activating enzyme (e.g. Wasp 
or ActA), actin polymerizes into a dense crosslinked filament 
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FIG. 3: Schematic view of the evolution of an actin gel (shades of 
red) around a bead (grey) as described in the text with (i) a sym- 
metric growth, (ii) the symmetry-breaking, and finally (iii) the comet 
formation. 



network, which extends several pm into the cell and where 
fast polymerizing barbed ends are oriented towards the mem- 
brane. Polymerization is restricted to a narrow zone close to 
the membrane since free barbed ends are rapidly blocked by 
capping proteins. Free pointed ends which are far away from 
the membrane are depolymerizing. These out of equilibrium 
growth processes driven by the irreversible hydrolysis of ATP 
lead to the extension of membrane protrusions or propel bac- 
teria forward. 



B. Biomimetic experiments 

1. General obsen'ations 

The minimal set of biochemical ingredients to induce actin 
driven motion have been identified about ten years ago llssll . 
At the same time actin driven motion has been successfully re- 
constituted in vitro by re plac ing the bacterium by mimetic ob- 
jects, e.g. beads i[Tlll4lll62ll . vesicles {23, 35, 56] or droplets 
a. 

In these experiments objects (hard, soft, fluid) coated with 
either ActA or Wasp/Scar proteins are added to a solution con- 
taining ATP actin, capping protein and a few well defined reg- 
ulatory proteins. With this design, actin polymerizes predom- 
inantly at the surface of the object, i.e. the internal interface, 
and depolymerizes at the interface between the network and 
the solution, i.e. the external interface. 

After an initial phase (i), where polymerization occurs sym- 
metrically around the object, the symmetry is broken and the 
actin cloud starts to grow asymmetrically (ii). In later stages 
an actin comet develops (iii), and the object starts to move 
with velocities up to 0.1 |jms~^. A schematic representation 
of the three phases of the actin cloud evolution is shown in Fig. 
[3] Interestingly the mode of movement depends on the surface 
parameters of the object. Depending on the conditions one can 
observe a continuous motion or a saltatory motion, where the 
object undergoes stop-and-go cycles isIlSlllllSSl], which is 
also reflected by variations in the actin density in the comet. 

The network grown around these biomimetic objects has 
elastic properties with a Young's modulus of lO'^-lO^'Pa ll20l 
[36ll and is often referred to as an actin gel. 

The actin filaments interact, at least transiently, with the 
activating enzyme bound to or adsorbed on the objects sur- 
face, e.g. the stress to detach an actin comet from a bead has 
been estimated to be of the order of ~ 100pN|jm~^ ll36ll . 
whereby the adsorbed activating enzyme stays on the bead. 



Interestingly, in the same experimental setup actin comets un- 
der compression appeared to be hollow. 

In regime of continuous motion the bead velocity is not 
affected by the viscosity of the medium (over five orders of 
magnitude). This raises the question of the dissipative force, 
which is obviously not the Stokes force on the object (about 
1 fN for a bead of 1 |jm radius moving with a velocity of 
Ifimmin^^ in a standard motility assay). In contrast, the 
stalling force for the growing actin tail is a few nN ll36ll48ll ). 
Several studies indicate that friction between the actin gel and 
the propelled object is the major source of dissipation lfl9ll6lll . 
This would support the hypothesis, that the saltatory bead or 
vesicle motion is the result of a stick-slip motion, where a 
certain critical force has to be overcome to rupture bead-gel 
bonds and to displace the object with respect to the gel ll4l [l9tl . 



2. Symmetry-breaking 

Biomimetic experiments are not only an effective tool to 
study the generation of motion by polymerization but they 
have also revealed a spontaneous symmetry -breaking instabil- 
ity in the growing actin network (transition form (i) to (ii) in 
Fig. O, whose consequences in vivo are not clear and whose 
nature is still a subject of debate in the literature. 

Even though, as often argued by biologists, a symmetry- 
breaking does not play a role in bacterial systems, since the 
ActA proteins are distributed asymmetrically around the bac- 
terium due to a previous cell division, the instability is a pow- 
erful tool to study the coupling mechanisms between actin 
polymerization and mechanical stresses. 

So far there seems to be a consensus in the literature, 
that mechanical stresses build up in the actin network due 
to growth, since new monomers are inserted at the internal 
curved interface and push older network layers away from 
the object. In a symmetric situation one expects therefore 
high tangential stresses at the external actin/solution interface 
and high normal stresses at the internal object/actin interface. 
Then the symmetry-breaking is driven by a release of elastic 
stresses in the actin gel, either by an asymmetric polymer- 
ization/depolymerization or by a fracture at the external gel 
interface, whereby the two mechanisms are difficult to distin- 
guish. 

There seem to exist subtle differences in the nature of the 
syrnmetry-breaking for hard ifTll l4ll l62ll and soft objects 
ll23i[35ll56ll . as was demonstrated in a more recent study ll4ll . 
In fluorescence labeling experiments it was shown that for soft 
objects (vesicles) the symmetry is broken at the internal gel 
interface, i.e. polymerization is considerably slowed down on 
one side of the vesicle, such that an actin tail develops at the 
opposite side. In contrast, for hard beads the symmetry is bro- 
ken at the external (depolymerizing) gel interface, whereby 
some authors suggest, that the filament network is ruptured 
due to an accumulation of tangential stresses ifssll . whereas 
theoretical models indicate that a local stress dependent de- 
polymerization is sufficient to induce a symmetry-breaking 
l28l [soil . Note however, that actin driven motion itself does 
not require the polymerization at curved surfaces, as has been 
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demonstrated theoretically lfT2ll . 

Furthermore, the mechanical properties of the biomimetic 
object has an effect on the mode selection of the instability. 
Whereas for soft objects only instabilities, which produce one 
single actin tail, have been reported, there are observations of 
higher order instabilities for hard spheres depending on the 
experimental conditions ifTill . Therefore, the boundary con- 
ditions at the internal gel interface seem to be crucial for the 
mode selection. 



A. Brownian ratchet models 

I. General concept 

As mentioned above, the idea of a Brownian ratchet model 
for polymerization forces in biological system was introduced 
by Peskin and coworkers ||44}l . This work explored the recti- 
fication of the Brownian motion of a particle by the intercala- 
tion of new monomers at the interspace between a filament tip 
and the obstacle, which gives an ideal ratchet velocity v of 



III. THEORETICAL APPROACHES 

Two types of theoretical concepts, how polymerization pro- 
cesses transform biochemical into mechanical energy, have 
been developing in parallel over the past 15 years. 

The first type of concept, called "Brownian ratchet models" 
was introduced by Peskin et al. lliill . Its major postulate is that 
polymerization processes (e.g. actin polymerization) are able 
to rectify the Brownian motion and can thus induce motion. 
Brownian ratchet models describe microscopically the poly- 
merization of actin filaments in the presence of an obstacle. 
Although very consuming in terms of computational efforts, 
Brownian ratchet models allow for the incorporation of a very 
detailed description of the kinetics of the polymerization ma- 
chinery. They provide ingenious tools to study the complex 
phenomena which have been observed in biomimetic experi- 
ments. 

The second type of concept, we shall refer to it as macro- 
scopic concept, does not focus so much on the dynamics of the 
single filament, but rather considers the actin filament network 
as a continuous elastic body under growth, where the growth 
dynamics is driven by a thermodynamic force, the chemical 
potential. These coarse-grained models emphasize the global 
stress distribution in the filament network and the nonlocal 
aspect of elasticity, while neglecting the details of the poly- 
merization process. Nevertheless, they are more suitable to 
describe and understand the symmetry-breaking in actin gels 
around spherical objects in terms of simple physical ingredi- 
ents. 

The ideal multiscale model would combine both concepts 
and use the global stress distribution as an input for the com- 
plex polymerization kinetics at the free interfaces. 

In the following section we will briefly review the two con- 
cepts. In the first subsection on Brownian ratchet we will 
mainly discuss Refs. lfl6ll2lll as they provide, from our point 
of view, the most advanced description for the polymeriza- 
tion dynamics close to the obstacle. In the second subsection 
on macroscopic models we will shortly outline and discuss 
the major drawbacks of Refs. jsotl and |j28|] and then give a 
perspective, how these problems can be solved using homog- 
enization techniques. 



where D denotes the diffusion coefficient of the obstacle and 
5 the size of a monomer. However, later on it was shown, 
that the actin driven motion was relatively independent of 
the obstacle size ||24i 1 and independent of the viscosity of the 
medium over several orders of magnitude ll6lll . i.e. inde- 
pendent of D. Therefore, the concept of Brownian ratchets 
was generalized to "elastic Brownian ratchets", where the tips 
of polymerizing filaments undergo fluctuations, which induce 
a propulsive force on the obstacle |[39l [59tl and to "tethered 
ratchets" ll40ll . which involves also the transient attachment of 
fluctuating filament tips to the obstacle. Typically in this de- 
scription the polymerization rate constant fcp is weighted by 
the load / using Kramers theory 

where fc™"^ is the rate constant at zero load, 5 represents the 
gap size to intercalate a monomer, fc^ and T denote the Boltz- 
mann constant and the absolute temperature, respectively. 

The general framework has been used to quantitatively 
model the steady motion of flat objects lfl2ll . lamellipodia 
and bacterial motion ll39il40ll . However, stochastic effects in 
the number of polarizing filaments were necessary to break 
the symmetry in a spherical gel around a bead homoge- 
neously covered by ActA or Wasp, whereas the global elastic 
stress distribution due to filament crosslinking was neglected 

Other attempts to model biomimetic motility quantitatively 
by explicitly modelling the filament dynamics in the actin tail 
succeeded in obtaining the crossover from a continuous to 
a hopping motion |^ . However, the velocity oscillations 
were on a time scale of ms with stepsizes of a few nm, as op- 
posed to the experimental oscillations on the scale of several 
min with stepsizes on the pm scale lisi [Tii [sillssll . 

More recently, the tethered ratchet model of Mogilner et al. 
ll40ll was submitted to a more rigorous treatment concerning 
the poly mer physics of the filament brush close to the obsta- 
cle iillll. Further away from the obstacle the cross-linked 
gel is advancing with a so-called grafting speed (force depen- 
dent and coupled to the brush length). In the following we 
will highlight the basic features of this model, since it quanti- 
tatively reproduced velocity oscillations. 
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FIG. 4: Schematic representation of the elastic ratchet model as ex- 
plained in the text. Filaments are anchored with one end into a cross- 
linked network forming an actin-gel and oriented with the other end 
against an obstacle. Adapted from Refs. fT6lf2lll . 

2. A quantitative model for velocity oscillation in actin-based 
motility 

As mentioned previously, the motion of a mutant form of 
the bacterium Listeria (with a mutation in the ActA protein) is 
oscillatory and shows remarkable temporal patterns lfl9l[3lll : 
the bacteria move very slowly during 30 to 100 s, jump for- 
ward during a few seconds and then slow down again abruptly. 
Such a periodic behavior, consisting of long intervals of a 
slowly changing dynamics that alternate with short periods of 
very fast transition, are found in several chemical and biolog- 
ical systems and are known as relaxation oscillations |30|]. 

Gholami et al. ll2Tll and Enculescu et al. rtT6ll have devel- 
oped a microscopic model based on the concept of tethered 
elastic ratchets for actin-based motility. Their model consists 
of a brush of growing actin filaments close to an object (the 
bacterium) and describes the generation of forces and conse- 
quently the propulsion of the object. 

Briefly, the model considers the case of fluctuating fila- 
ments close to an obstacle. Filaments may attach to the obsta- 
cle with a rate constant ka and detach from the obstacle with 
a force dependent rate constant kd, resulting in two distribu- 
tions of populations of filaments, i.e. of attached and detached 
Ua and Ud, respectively. Opposite to the obstacle the filament 
ends are anchored in a cross-linked network, the actin gel. De- 
tached filaments polymerize with a load dependent velocity 
Vp. The distance between the grafting point, i.e. the interface 
between the network and the polymer brush, and the obstacle 
is denoted by ^. The filaments are characterized by their free 
contour lengths I. One of the crucial ingredients of the model 
is, that the grafting point is advancing in the direction of the 
obstacle with a so called grafting velocity Vg, which depends 
on the free contour lengths of the polymers by 

^'ffW =<""tanh (//[), (4) 

where I denotes a characteristic width of boundary between 
the crosslinked network and the filament brush. 

Attached and detached filaments exert entropic forces ll22ll 
on the obstacle, Fa{l, and Fd{l^ C)' respectively, which lead 
to the propulsion of the object with an effective friction coef- 
ficient 



The load dependence of the kinetic rate constants is again 
included using a Kramers type expression [see Eq. (O], i.e. 

kd{Fa)=k''dGM-5Fa/kBT) (5) 

for the detachment and 

t'p(0=«r"exp(-<5F<i/fcBr) (6) 

for the polymerization speed. 

The full evolution of the length distributions of the two 
filament populations na{l,t) and nd{l,t) is described by 
advection-reaction equations. However, it is shown that the 
two distributions contract rapidly on the scale of 10~^ s into 
monodisperse distributions Na = na{t)5{l — la{t)) and Nd ~ 
nd{t)5{l — ld{t)) localized at la and Id for attached and de- 
tached filaments, respectively. Therefore, the dynamics of the 
system can be simplified to four ordinary differential equa- 
tions for the evolution of the free contour lengths /^(t) and 
ld{t), respectively, the number of attached filaments na{t) 
(the total number of filaments N is constant and therefore 
nd{t) = N — na{t)) and the distance ^(t) between the ob- 
stacle and the grafting point. 

The solution behavior of this system of equations has 
been analyzed numerically depending on the maximal graft- 
ing speed and the rate constant for the attachment of 
filaments ka- Fortunately, most other model parameters are 
known experimentally. The model displays two different dy- 
namical regimes: steady and oscillatory motion, whereby the 
oscillatory regime is robust against changes in the parameters. 
The oscillations occur on the time scale of min and produce 
jumps of the obstacle displacement in the pm range and re- 
semble very much relaxation oscillations. 

A deeper analysis of the solutions suggests, that the oscilla- 
tions arise from a so called push-pull mechanism, i.e. a com- 
petition between pulling and pushing forces acting on the ob- 
stacle. In this mechanism, a long pull phase, where most of 
the filaments are attached to the obstacle and polymerization 
stalls, alternates with a short push phase, where most filaments 
are detached and polymerize rapidly. In the pull phase the 
grafting velocity Vg is higher than the polymerization speed 
Vp, and the magnitude of the forces on the obstacle increases 
due to their dependence on ^, la, and Id- At a certain point, 
in the pushing phase, the pushing forces outweigh the pulling 
forces and cause an avalanche like detachment of filaments 
and the obstacle "hops" forward, lowering the load on the de- 
tached filaments, which start to polymerize rapidly and thus 
Vg < Vp. Meanwhile free filaments start to attach to the ob- 
stacle and increase the pulling force, i.e. the obstacle slows 
down. Free filaments start to buckle, the polymerization stalls 
and the cycle reenters the pulling phase. 

Obviously, this complex cycle arises from the subtle inter- 
play between pushing and pulling forces on the one hand and 
the grafting and polymerization velocities on the other hand. 
It would be an interesting task, to explore the parameter space 
a little bit further and identify the absolutely necessary ingre- 
dients to find oscillations in the push-pull mechanism. 

Besides steady and oscillatory motion the model also yields 
bistable and excitable behavior, which might lead to the rein- 
terpretation of previous experiments and is reminiscent of the 
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behavior caused by nonlinear friction in a variety of systems 
with a complex surface chemistry IstIi . Another aspect for 
future work is to couple the microscopic dynamics of the fila- 
ment brush to the bulk mechanics of the cross-linked gel. 



B. Macroscopic models 

While microscopic models give a detailed description of 
the polymerization and cross-linking dynamics, macroscopic 
models are more concerned about the global stress distribution 
in the gel, but adopt a more general formulation for the inter- 
face dynamics. In the following we will shortly introduce and 
discuss the essence of three simple models, i.e. by Lee et al. 
iH, Sekimoto et al. ISO:], and John et al. describing the 
symmetry-breaking in an actin gel around a solid bead. We 
will conclude this section with the discussion of a more ad- 
vanced mechanical model using homogenization techniques 
proposed by Caillerie et al. idT 



1. A phenomenological model of symmetry-breaking 

As an introduction we will present the problem of 
symmetry-breaking from a purely phenomenological point of 
view as proposed by Lee et al. ll32ll . If a bead moves, it is 
natural to assume that a force is applied on the bead, proba- 
bly due to deformed filaments that releases stress on the bead. 
Let g(r, t) denote the force per unit area in the normal direc- 
tion that is exerted on the bead, which could be e.g. a function 
of the local actin concentration. The total force on the bead 
is given by F = J ^(r, t)ndA, with n being the unit normal 
vector on the bead. The velocity of the bead is related to the 
total force via a linear relation 



v^CF^^y" 9{r,t)ndA, 



(7) 



where ^ is a dissipative coefficient, taken to be scalar for sim- 
phcity (it is necessarily so for a sphere in a Newtonian fluid). 

It is then assumed that the rate of change of g{v,t) is a 
local function of the bead velocity v, and that there is a feed 
back of the motion on the force g: faster motion is associated 
with a decrease of polymerization in the front and an increase 
of polymerization at the rear, and hence has an impact on g. 
Under the assumption of analyticity, the evolution of g can be 
written as 



dtg 



eg 



5gv.n , 



(8) 



with c > in order to ensure stability of the homogeneous 
stationary state g = 0. The signs of a and b are left arbitrary 
for the moment. Eqs. ([8]l and O constitute a complete set that 
can be solved numerically or analytically using a perturbation 
ansatz as will be outlined in the following. 

It is a simple matter to see that the set [Eqs. (O and (|7]i] 
admits the fixed point g = v = 0. By superposing small per- 
turbations on this solution, reporting into the above set, and 
expanding up to linear order in the perturbations, one finds 



that the fixed point is unstable for a > 3/^ = Oc and sta- 
ble otherwise. If motion takes place then this means that a 
symmetry-breaking has occurred, and thus the local force g 
has lost the spherical symmetry. More precisely by expanding 
g in spherical harmonics, g = J2em 9^mYem{9,(t)), and re- 
porting into ^ by assuming a direction of motion, say along 
oz, one finds to leading order 



dtg in 



1 



-Sn 



gtm + 



(9) 



whereby the ". . ." refer to nonlinear terms. Integration of 
other harmonics than the mode I ^ I'm^ vanishes exactly 
due to symmetry. Beyond the symmetry-breaking bifurcation 
for a > 3/^ the perturbations grow exponentially in time, and 
nonlinear terms are needed. Since only the mode I — 1 \s 
excited, the other modes are treated as adiabatically enslaved 
to it (at least in the vicinity of the threshold). The following 
analysis consists in expanding the solution for higher order 
terms and expressing the amplitudes of the higher harmonics 
in terms of the amplitude of the first order mode. Once g has 
been replaced by v from (|7]) and inserted into ^ a closed 
equation for v (which is simpler to assess experimentally than 
g) is obtained Il32ll 



9tv = ev + —7 



27 

5^2 



3 



(10) 

e = (1 — a/ttc) is a small parameter (expressing the fact that 
we focus on the instability threshold), and we have used the 
convention v'^ = v^w, and so on. The expression of u in terms 
of the coefficients a, 5, c, and ^ is not shown here. However, 
it is reported in |32] that u is negative in the considered pa- 
rameter space. The amplitude equation ( fTOl i serves to discuss 
the phenomenology of the motion, that we will briefly sum- 
marize here and which is shown schematically for the h — a— 
parameter plane in Fig. |5j 




FIG. 5: Phase diagram arising from the analysis of Eq. dlOb . The 
transition between the stationary (white region) and moving (light 
shaded region) state occurs at a = ac. For hd < b < bc2 the 
bifurcation is supercritical. In the dark shaded region there coexist 
two solutions v = and v = const 7^ which are linearly stable. 
Adapted from Ref (H. 

If the coefficient of the cubic term is negative, then a contin- 
uous transition from the symmetric (motionless) to the asym- 
metric (moving) state occurs at e = 0, i.e. the bifurcation is 
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FIG. 6: Schematic view of a bead surrounded by an elastic gel, show- 
ing several definitions as explained in the text. Shown is a symmetric 
gel (left) and an asymmetric gel (right). In the latter case the ten- 
sile stress is rather generally denoted by arr with r designating the 
tangent direction. 

supercritical (an analogue of a second order transition). This 
happens for a certain range of the product £fi and by fixing c 
to unity for definiteness. For a certain range of the product 
^6 the cubic coefficient changes sign, a signature of a sub- 
critical bifurcation (an analogue of a first order transition). 
This is known to lead to multistability: in a certain parame- 
ter range for e < there coexist two solutions: v = and 
V ~ const, both solutions being locally stable, i.e. stable with 
respect to small perturbations. One may think that the system 
may switch back and forth between the two solutions. 

Several remarks are in order (i) The above discussion 
is fully phenomenological, and actually Eq. ( fTOl i could have 
been written directly. However, the derivation given in ll32ll 
has a certain merit with regard to the description of the gen- 
eral laws behind the motion, and the feedback mechanism, 
(ii) The coefficients are certainly complicated functions of ex- 
perimental parameters, and a microscopic model is needed in 
order to relate phenomenological and experimental parame- 
ters, (iii) The origin of the forces acting on the bead and the 
subsequent generation of motion in this model is not clear We 
will come back to this vital question in the final conclusions. 

2. The role of tensile stress during the symmetry-breaking in actin 

gels 

A first physical macroscopic model including elasticity has 
been put forward by Sekimoto et al. ifsoll in order to explain 
the birth of the symmetry-breaking of an initially symmetric 
gel layer growing on a spherical or cylindrical object. The 
crux of their analysis is that the gel that has been formed is 
being continuously pushed outwards due to the arrival of new 
monomers at the bead. This is supposed to lead to large lateral 
stresses on the outer gel interface. A schematic representation 
of the model is shown in Fig.|6] In the following we outline the 
model in a cylindrical geometry. If age designates the tangen- 
tial stress, and if the gel layer is axisymmetric, it is postulated 
that 

r — rn 

<Jgg = E ^, (11) 

ro 

where E is the Young modulus. At r = (the bead ra- 
dius) the tensile stress vanishes, expressing the assumption 



that the gel layers are added unstretched at the bead surface. 
If h denotes the total gel thickness, one can write for the ten- 
sile stress at the external gel surface <7gg\r=r„+h 
The external surface of the gel is not subject to a force, so 
that both the tangential and normal forces must vanish, i.e. 
arg\r=ra+h = cTjr |r=ro+/i = 0. The tangential forcc is also 
zero at the bead surface, arg\r=ro- For a symmetric gel the 
shear stress vanishes everywhere in the gel. If a non zero ten- 
sile stress exists [Eq. (fTTTil. the mechanical equilibrium in the 
bulk, i.e. 

div{a) = (12) 

provides the following relation between the tensile and radial 
stress distribution, agg and arr. 

The gel may grow because of a gain in polymerization energy 
at the bead surface. This growth takes place at a certain price: 
the higher the thickness of the gel is, the larger is the stored 
elastic energy. Therefore, one expects the growth to stop at 
a certain equilibrium thickness Hq. The following growth ki- 
netic relation has been suggested ifsoll 

dth = kpe''^''-''\-=-« - fc^e=^'^«<'l-='-o+'. (14) 

where kp, kd, Cp, Cd are positive constants. The first term, 
which is positive, accounts for the polymerization at the bead, 
while the second one, which is negative, refers to depolymer- 
ization at the external surface. Note that on the one hand 
(Jrr\r=ro = —Eh'^/{2r'^) < 0, and thus stress penalizes 
polymerization and acts against gain in chemical bonds at the 
bead. On the other hand <Jgg\r=ro+h > and this causes 
the depolymerization to increase with the gel thickness, since 
(Jee\r=ro+h increases with h. Setting dth = provides us 
with the steady state thickness ho as a function of other pa- 
rameters. One straightforwardly finds 

+ J_i^:^) = o. (15) 

Cd ro ro CdE kd 

A steady solution exists as long as kp/ kd > 1- It can easily be 
checked that the steady solution hg is stable with respect to a 
homogeneous, i.e. axisymmetric, increase or decrease of ho. 
One further sees that ho is a linear function of ro, since only 
the ratio ho /ro enters the above equation. This seems to be in 
good agreement with experimental observations ll4lll . 

Let us now discuss the linear stability analysis of the steady 
solution with respect to perturbations that break the circu- 
lar symmetry. In a cylindrical geometry the eigenmodes are 
~ e""^o (qj- cos(m0o) in real variables) where 9o is the angu- 
lar variable shown in Fig.|6] and m is an integer Because the 
equations are autonomous with respect to time, all the eigen- 
modes can be written as ~ e^*, where (3 is the amplification 
or attenuation rate of the perturbation that must be determined 
from the model equations. We set for the perturbed thickness 

/i(0o) = /io[l + em(t)cos(m0o)]. (16) 
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e„i is a small, time-dependent quantity that justifies the linear 
stability analysis. The stress field inside the gel will thus be 
parametrized by the function h{6o). The linear analysis now 
consists in (i) solving the stress field in the gel with the appro- 
priate boundary conditions, and (ii) using the kinetic relation 
(fl4l i to obtain the dispersion relation /3 = /(m, p), where p is 
an abbreviation for all the physical and geometrical parame- 
ters that enter into the equation. 

In order to solve for the elastic field, one needs tospec- 
ify a constitutive law. Unlike in the two models ifiol l28ll we 
will discuss in the following sections, where a constitutive law 
is used by evoking basic continuum mechanics concepts, the 
Sekimoto et al. model ifsoll is based on an extension of the 
postulate represented by Eq. ( fTTT i to a modulated thickness 
/i(6'o). The basic ingredient of their analysis is to introduce 
an unknown function 0{6o) such that a material point of the 
gel which is originally located at Oq in an undeformed refer- 
ence state is moved to a new position 9(9q) upon deforma- 
tion. Then the elongation ratio (r — ro)/ro is replaced by 
{rd6{9Q) — rodOo) / {rodOo) so that the tensile stress takes the 
form 



= E 



.) 



ro d9o 



- 1 



(17) 



By using the equilibrium balance condition Eq. ( fT2] i and ne- 
glecting the shear stress (for a discussion of this assumption 
see Ref. llsoll ) one can again find a relation between arr and 
aee- By defining 



one deduces that 



pro+h 

T = (Jggdr 

■I ro 



T 



(18) 



(19) 



since only this quantity enters the kinetic relation ( fT4l i that is 
needed for the derivation of the dispersion relation. Recall 
also, that (Trr\r=ro+h = to linear order in the perturbation 
Em- Upon using ( fTTj i the integral ( fTSl l leads to 



T = E 



h{9o) 



h{9o)^\ d9{9o 



2rn 



d9„ 



- h{9o] 



(20) 



Additionally, Eq. ( fT2] ) yields, under the assumption of zero 
shear stress, deT = 0, that is T is independent of 9. This 
condition provides a relation between 9{9o) and h{9o). Upon 
using J^^ ^Iw^'^^o = 27r, one can express T as a function 

of an integral J^^ F{h{9o)) d9o where only h{9o) enters (F 
is given by Eq. (D5) of Ref. llsoll ). Plugging this relation 
into ( [2Q] i provides a relation between d9/d9o and h{9Q), and 
substituting hhy Eq. ( fT6l ) (Tgg\r=ro+h and crrr\r=ro can be de- 
duced to first order in ■ After some algebraic manipulations 
the dispersion relation is obtained 



/3 



ro 



(21) 



solution 
p=0 




FIG. 7: Schematic view of a bead surrounded by an elastic gel, show- 
ing also several definitions of Ref. 12^ as explained in the text. 



where 



nr. 



CdEe 



CdEha 



ho 



ho 



(22) 



and where we have set /lo/^o = ^o- P is positive, meaning 
that the perturbation grows exponentially with time: the sym- 
metric gel layer is thus unstable. Surprisingly, the dispersion 
relation does not depend on m (since ^lm does not). Conse- 
quently, all wavenumbers have the same growth rate. Thus, 
the linear stability analysis does not select a typical mode m 
(like the fastest growing mode) for the instability. 



3. A nonlinear study on symmetry-breaking in actin gels 

Unlike the previous model where a tensile stress distribu- 
tion is a priori postulated, the idea of the model proposed by 
John et al. ll28ll is to treat the actin gel as an elastic continuum 
in the framework of a linear theory, and to formulate a simple 
kinetic relation expressing growth (or polymerization), differ- 
ent from Eq. (T4[ . 

The model considers a bead (radius ro) surrounded by a 
growing elastic actin gel (radius ro + h) as shown in Fig. |7] 
The gel is stressed by a small molecular displacement L in 
normal direction at the bead/gel interface, i.e. 



(23) 



where Uj. denotes the radial component of the displacement. 
This choice is motivated by the microscopic picture, that for 
the addition of monomers, enzymes facilitate a molecular dis- 
placement L at the bead/gel interface. This displacement is 
the source of stress. The bead/gel interface as well as the ex- 
ternal gel surface is shear free and the normal stress at the 
external surface is set to zero^. No condition is imposed on 
the normal stress component at the bead, since there a dis- 
placement is imposed instead. The boundary conditions on 
the stress are thus 



rti(JijTj\r—rQ,r—rQ + h — 
rii(^ijrij\r—ro+h — : 



(24) 



^ Actually it can be set to —p, where p denotes the liquid pressure, but this 
contribution is quite small. 
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where rii and are the ith component of the unit normal and 
tangent vector of the surface under consideration (bead or ex- 
ternal gel surface). 

The stress distribution in the gel is obtained then by solving 
the Lame equation for the displacement field 



1 



1 - 2s 



V(V.u) = 0, 



(25) 



where s is the Poisson ratio. The stress is related to the dis- 
placement u by Hooke's law 



(26) 



where e.y = {diUj + djUi)/2 is the strain tensor, and A and 
1/ are the Lame coefficients which are related to the Young 
modulus E and s (for an isotropic material there are only two 
independent elastic parameters). 

Eqs. ( |23] | - dZST l represent a complete set that allows to deter- 
mine, the stress and displacement fields in the gel. Note that, 
despite the fact that the bulk equations are linear, the problem 
acquires a nonlinear character via the geometry of the exter- 
nal gel boundary. Indeed if we fix a certain arbitrary geome- 
try h{9), then the stress and displacement will be a nonlinear 
function of h. The calculation can be handled analytically 
for a symmetric gel as well as in the linear stability analysis 
ll28ll . Beyond a linear analysis, a numerical study has been 
performed and will be briefly discussed below. 

Once the mechanical problem is solved, one needs to com- 
pute the cost in elastic energy per unit mass for inserting a 
monomer on the bead/gel interface, i.e. the chemical poten- 
tial difference. For the sake of simplicity we focus on the 
case of an inert external gel surface, i.e. neither polymeriza- 
tion nor depolymerization takes place at the external interface. 
We assume that insertion of a monomer at the bead results in 
a displacement of the external gel surface in the radial direc- 
tion proportional to the chemical potential change at the bead. 
We will critically assess this assumption in the next section on 
homogenization models. 

One may then write a kinetic relation of the shape evolution 
of the gel envelope 



dth = -MAfi , 



(27) 



where M denotes a mobility and A/i the difference in the 
chemical potential between a volume element in the gel and in 
solution at the internal interface. Here we assume that the mo- 
bility is associated with the polymerization/depolymerization 
kinetics, which constitutes the prevailing dissipation mecha- 
nism. The chemical potential is composed of a contribution 
due to the gain in polymerization (denoted as A/ip < 0) and 
an elastic part ll28ll 



A/i = A/ip + vuijUi 



A 



FIG. 8: The dispersion relation as a function of the mode I obtained 
from the linear stability analysis of the John et al. model (2^ . The 
fastest growing mode corresponds to ^ = 1. 



stress combination. Eq. ( fT4l i contains only a linear form and 
no quadratic forms as in Eq. ( |28] l. Here, the quadratic form is 
essential for the mode selection leading to a comet formation, 
as will be described below. 

The stress problem Eqs. (|23]l-(|25Tl can easily be solved an- 
alytically for a spherical geometry (axisymmetric growth). 
Upon setting dth = {) in Eq. dZT]) one finds a steady solution 
with a gel thickness Hq obeying luSll 



^Ea - (1 - 2s)A/ip\ 
'2£;a+ (l + s)A/(pj 



1/3 



- 1 



(29) 



'■kk 



,{l+Ukk) (28) 



where a = i/f'o- This solution exists for 2Ea/{l + s) > 
— A/Xp! elasticity acts against monomer addition, so that 
the gel stops growing at that thickness. In the opposite 
limit growth continues without bound. Both situations have 
been observed experimentally ifisll . however in the latter case 
growth stopped due to monomer depletion at the bead/gel in- 
terface. 

The linear stability analysis of the symmetric case can be 
performed analytically (by decomposing the stress and the 
shape evolution into spherical harmonics Yi,n). The disper- 
sion relation (3{t) is presented in Fig. [8] The basic result is 
that a symmetric shape is unstable against symmetry breaking 
ll28ll . Interestingly the mode which corresponds to a transla- 
tion of the external surface with respect to the bead is the most 
unstable*. This translation motion is similar to that shown in 
Fig.|9](b). For this instability the quadratic terms in Eq. ( |28l ) 
play a crucial role, since considering only the linear terms 
leads to a stable symmetric solution, with a zero growth rate 
for the translational mode. 

In order to ascertain the subsequent evolution of the exter- 
nal boundary (i.e. in the fully nonlinear regime), a full numer- 
ical analysis has been performed ll28ll . The set of mechanical 
equations ( |23] l - ( |25l ), and the growth kinetics ( |27l ) have been 
cast into a phase-field approach, which has now become a fre- 
quent method to treat free moving boundary problems. For the 
details of the phase-field formulation and their numerical im- 
plementation we would like to refer the reader to the original 



Note that (|27] | differs from ( fT4l i not only by the presence of the 
exponential function (which can be linearized since the stress 
energy is always small in comparison to the thermal excitation 
energy; hidden in the constant Cp), but most importantly by the 



* We should not confuse the mode t. = 1 with the usual global translation, 
which is a neutral mode. Here only the external gel surface moves while 
the bead is fixed, so that the mode ^ = 1 is a physical one. 
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paper ll28ll . Here we shall only report the basic results 

The numerical study in two dimensions with plane strain 
shows that for axisymmetric initial conditions with small am- 
plitude perturbations, the symmetry is broken for the mode 
m = 1 (identical to the mode £ = 1 in three dimensions), 
which corresponds to a translation of the gel layer with respect 
to the bead. Where this symmetry-breaking occurs, depends 
only on the initial conditions. The instability then evolves fur- 
ther into an actin comet, reminiscent of the comet developed 
by Listeria monocytogenes. Fig. |9] shows a typical result of 
a numerical simulation. The comet formation seems here to 





□ 





FIG. 9: Symmetry-breaking of a circular gel. Shown is the evolution 
of the gel thickness, starting from a homogeneous thin gel (a) with 
random small amplitude perturbations, (b) shows the initial symme- 
try breaking, while (c) shows the subsequent evolution of the shape 
into a comet in the far nonlinear regime (adapted from Ref. ll28ll ). 

be the generic growth mode. This finding points to the fact 
that the comet formation is probably a quite robust feature; it 
results from simple physical prototypes. 

The physical picture of the symmetry breaking may be un- 
derstood as follows. Let us start with a symmetric layer as 
in Fig. |9] (a). Suppose that due to some natural (inevitable) 
fluctuation, the gel layer becomes asymmetric, as in Fig. |9] 
(b). The stress is due to addition of new monomers at the 
bead/gel interface. On the side where the gel thickness is 
small the stress field is stronger than on the other side, since 
the gel feels more the "outwards pushing" of new monomers 
inserted at the bead. Because of the increase of the stress 
(and strain) in the thinner layer, polymerization becomes un- 
favorable there. New monomers will preferentially be inserted 
on the side where the thickness is larger The appearance of 
modes larger than i ~ 1 would create several thin and thick 
regions, which are likely unfavorable. It seems thus that the 
mode £ = 1 is optimal for the insertion of monomers at the 
bead/gel interface. 

By considering also the stress dependent depolymeriza- 
tion at the external gel interface, one finds another instability, 
whereby the location of the most unstable mode is determined 
by surface tension. This result is in agreement with the occa- 
sional experimental observation of higher order modes 13, 
which will be discussed next within a more elaborate homog- 
enization model. 



4. Homogenization models 

So far, continuous models have been suffering from the in- 
adequate or insufficient description of the mechanical aspect 
of the actin gel. The greatest difficulty in the description of the 
mechanical equilibrium of the growing actin network arises 
from the fact, that the growth history determines the network 
structure and therefore also the stress distribution. Conse- 
quently, a realistic model would have to include the informa- 
tion on how the network evolved. 

A second problem, which is most prominent in the de- 
scription of the symmetry-breaking around solid objects, is 
the coupling of the growth process between the internal and 
the external interface. As an example, consider the problem 
of actin growth around a bead functionalized with ActA or 
Wasp/Scar; Typically growth takes place at the barbed end 
of the actin filaments, which points toward the bead/gel in- 
terface. This growth process pushes older gel layers further 
away from that interface. Experimentally, the growth process 
is observed by an increase in the gel thickness. However, it 
is not clear how the insertion of mass at a (fixed) solid/gel in- 
terface translates into the displacement of the (free) gel/liquid 
interface. A realistic elastic theory would have to account for 
this coupling problem in a rigorous way. 

One observation, which might help to solve the above men- 
tioned problems at least partly, is that the actin network forms 
more or less regular structures, which are not perfectly peri- 
odic, but could be considered in a first approximation as "al- 
most periodic". The actin gel can then be regarded abstractly 
as a network of elastic filaments connected by nodes with a 
certain periodicity. This network is completely defined by the 
positions of the nodes and their connectivities. In the network 
structure, the size of each elementary cell, e.g. the distance 
between two Arp2/3 crosslinks (~ several tens nm) , is small 
compared to the total size of the structure (~1 pm), which in- 
troduces a small parameter 77 into the problem, which is the 
ratio of the length of the unit cell and the total network size. 
In the following we will shortly outline the basic idea of an 
actin homogenization model in two dimensions ifioll . 

We consider a planar network of stiff elastic bars around 
a solid cylinder (radius ro) with the topology shown in Fig. 
[TOl The bars are connected to each other via nodes. The actin 
filaments are assumed to be linked to the cylinder at Nt sites 
evenly located on its surface at a distance p ~ ryrp between 
two close sites, i.e. at an angular distance 7] ~ j^. The actin 
gel is made of Nn layers of bars in the radial direction. As 
the growth of the gel is due to the polymerization of actin 
monomers at the surface of the cylinder, each layer is assumed 
to be made up of the same number of nodes. So the nodes 
of the gel can be numbered by two integers (i^^, i^^) with 
numbering the radial layers and 1^^ the position of the node in 
each layer, respectively. 

It is assumed that the discrete net is made up of a large 
number of bars meaning that Nt and iV„ are very large and 
of the same order To be more precise, the parameter rj is as- 
sumed to be very small and the number Nn of layers is given 
by Nn ~ — with a being of order 1 with respect to ?y. Using 
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FIG. 10: Sketch of a small part of the filament network with some 
examples for the node numbering (u^,u^). The "curves" along which 
i/^=const. and i/^=const. are shown in blue and red, respectively. 



this notation, a node of the gel can be labeled by (/i 



with /i*'' = rju^ The coordinates [/x^'', /x^'') take values in 
u! =]0, q;[x]0, 27r[ and are meant to become the set of La- 
grangian curvilinear coordinates of the equivalent continuous 
medium. 

The upscaling of the net to a continuous medium consists in 
determining the equivalent stresses from the bar tensions, the 
equations of equilibrium (or motion) satisfied by these stresses 
and an equivalent constitutive equation ensuing from the prop- 
erties of the bars. This can be carried out by using an asymp- 
totic expansion (for an introduction see ||^[33,ll3l)- Here, as 
the network structure is simple, a more heuristic presentation 
can be used. The basic idea of the homogenization process 
is that, for most of the motions of the network, the positions 
of its nodes (z^^, i^^) can be approximated by a continuous 

function ip (i/^, i^^) where /i*'' = 771/*. The purpose then is to 
determine the equations governing this deformation function. 
The equivalent Cauchy stress tensor is given by the classical 
relation due to Cauchy ll9il34ll 



1 ^ 

a = -yN' 



e ' 



(30) 



with B'^ being the "bar vector" linking the two extremities of 
the bar b [b ~ 1, 2, 3) of the elementary cell (shown in Fig. 



of the network in a deformed state, e 



being the corresponding unit vector, 



we'll 

N'' the tension in the 



bar and g 



being the surface of the elementary 

cell. The constitutive equation of the equivalent continuous 
medium follows from the constitutive equations of the bars 
which, for the sake of simplicity, are assumed to be 



N" 



F - r 



(31) 




FIG. 11: Sketch of an elementary cell, showing the node numbering 
{v^ , v^) and the elementary bar vectors (solid arrows). All other 
bars are shown as dashed lines. The "curves" along which !/^=const. 
and i/^=const. are shown in blue and red, respectively. 



In order for a symmetrical equilibrium configuration of the 
gel to be possible, the constitutive equations of the bars 1 and 
3 should be identical, i.e. ~ and l]-^ = Zj^. 

Following the homogenization assumption stating that the 
position of the node (i^^, i/^) is ip (^i'^, v^^ with /x"* = 'i]v\ a 
simple Taylor expansion yields: 



B^ = 
B^ = 
B' = 



rj ydiip + dl^tj} 



(32) 



where df = Since they are associated with a quite sim- 
ple numbering system for the nodes, the variables and /i^ 
arise naturally as Lagrangian variables of the equivalent con- 
tinuous medium through the homogenization process. How- 
ever, they are not the most convenient variables to study sym- 
metrical equilibrium configurations. Therefore we have intro- 
duced the variables and defined by 



Ai=^iandA2 = ^+/x2, 



Setting 



one finds 

B^^rjgi, ^2 ==7752 and B^ 



v -91 + 7752 



(33) 



(34) 



(35) 



with gi = and 92 = ^2 V'- 

Carrying these relations into ( [30l l yields 



with /J^is the length of the bar 6 at rest and serves 

as a parameter in the constitutive equation. 



1 



^S'® g. 



(36) 
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with 51 = 7] (iV^e 1 - N^e^) and S"^ 
7? (iiVie 1 + N'^e^ + iiV^e^). 

As the only forces acting on the gel are applied on its 
boundaries, the equilibrium of the continuous medium reads 
classically 

diver = 0. (37) 

Using the virtual power formulation of that equation and the 
change of variables 

it can be proven that the equilibrium equation reads 

2 

J2 diS' = . (38) 

i=l 

a. Coupling between growth and mechanics To study 
the growth of such an homogenized network, one can stay 
within the picture of a mechanical equilibrium of the actin gel 
on the time scale of the growth process. From the homoge- 
nized elastic equations one can derive the elastic contribution 
A^e to the chemical potential for the addition or subtrac- 
tion of nodes at the gel interfaces starting from the free elastic 
energy Ff. in the network, i.e. 

^ •'^ 6=1.2,3 

where is the elastic energy associated with the extension or 
contraction of each of the filaments. The elastic chemical po- 
tential is then given by the variation of the elastic energy with 
respect to the size and shape of the network A/ig = 5Fe/5Vl 
by respecting the boundary conditions. 

We assume now that the chemical potential contains also 
a contribution from the chemical process of polymerization 
A^c, where A/ic < at the internal interface and A/^c > 
at the external interface. This assumption accounts for the 
polar treadmilling behavior of the actin polymerization, i.e. 
polymerization occurs at the internal interface and depolymer- 
ization at the external interface. Furthermore, the chemical 
potential contains a contribution from interfacial energy, i.e. 
A/^s = — 7K, with 7 being the surface tension coefficient and 
K being the curvature of the interface. This leads to the fol- 
lowing expression for the normal velocities of the two free 
interfaces in the Lagrangian coordinates (A^, A^) 

v„ = -r^M' {Apl + Apl + Ai^i'^) = -ijM'Ap' , (40) 

with 2 = for the internal and i = 1 for the external interface. 

b. Homogeneous gel growth and linear stability analysis 
First one may consider the symmetric problem, i.e. the growth 
of a gel with homogeneous thickness a = rjNn which has an 
axisymmetric solution ~ (^r(Ai)e'r(A2). The equiUbrium 
equation in this case then reduces to 

= 28^ [N^d^^r] - QiVi + , (41) 
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FIG. 12: Filament lengths (top), radial and tangential tensions (bot- 
tom) depending on the positions in the network in mechanical equi- 
librium. Parameters are /k^ = 1 and l]-,-^ = = p. 




a 

FIG. 13: Dependence of the elastic part of the chemical potential 
A/ie on the number of radial filament layers a. Remaining parame- 
ters are — l^-,^ — p and /k^ — 1. 

where we have introduced 

~h N'' 

N' = il-jr- (42) 

In this situation the filaments with 6 = 2 are oriented in a 
tangential direction. Eg. (HTli can be solved numerically using 
continuation methods llTsIl . 

Fig. [T2] shows a solution of Eq. W\\ for a given network 
thickness a. Naturally, is extending as one moves away 
from the bead, whereas is first shortening and then extend- 
ing to reach its equilibrium length at the outer gel surface. 
Consequently, the gel is under radial compression and under 
tangential extension far away from the bead surface. However, 
for regions close to the bead surface the gel is under tangential 
compression. 

Fig.[T3]shows the dependence of the chemical potential de- 
pending on the number of radial filament layers a. Assuming 
that new filaments are inserted in the same stressed state as 
the already present material at the two interfaces, the elas- 
tic chemical potential is identical at the two interfaces for a 
homogeneous gel (A/Xg = A/i^ = Afig)- With increasing 
network size, i.e. increasing a, A/ig increases in a strongly 
nonlinear fashion. Note, that for higher values of a the ho- 
mogenization approach breaks down and the gel is starting to 
"fold back". Beyond this point, at a « 1.25 in Fig.[T3] no 
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4 6 

mode m 



FIG. 14: Dispersion relation with Dirichlet boundary conditions at 
the internal interface. Shown is only the larger of two eigenvalues /3i 
(the second one being always negative) depending on the wavenum- 
ber m for various values of the surface tension as indicated in the 



legend. Parameters are a = 
scale is r = ro/(Mk^p^). 



: p. 



The time 



curvature break the symmetry faster, than larger beads 15811 . 
which is in agreement with our model, where the time scale 
of symmetry-breaking increases linearly with the radius of 
the bead vq. Note also that, assuming a constant polymer- 
ization potential, but changing the radius vq by keeping all 
other parameters constant, leads to a linear relation between 
the gel thickness, i.e. (pr{a) — ipr{0), and tq in steady state. 
These two results have also been obtained in simpler models 
based on scaling arguments ll4lll48i Issll and hold for the case 
that monomer diffusion is fast enough to avoid depletion of 
monomers due to polymerization at the internal bead gel in- 
terface. Another interesting point is the type of instability one 
might observe, i.e. an undulating vs. peristaltic instability. 
For small modes m < 4 one finds an undulating instability, 
i.e. the perturbations at the external and internal interface are 
in phase, whereas for higher order modes on should observe 
a peristaltic instability, where the two perturbations have the 
opposite phase (data not shown). 



physical meaningful solutions exist. 

We consider now a filament network which is allowed to 
grow symmetrically with identical mobilities (M" = A/^ = 
M) at the two interfaces following the two dynamic equations 

dta° = T^M + A/ie) = -Vp + V, (43) 
dta^ = -rjM (A/iJ + A^e) = -Vd - Ve (44) 

whereby the positions of the internal and external interface 
are denoted by a° and a^, and where we have introduced 
the polymerization speed Vp = —rjMAfj!^ > 0, the de- 
polymerization speed Vd = yyMA^J and an "elastic speed" 
f e = rjM A/if,. The steady state for the gel thickness is given 
by (9ta = dta^ — dta^ = 0. Defining now a mean velocity 
V = {vp + Vd)/2 and a velocity difference Af = {vp — Vd)/2, 
one obtains that in the steady state Av = Ve and dta^ = 
dta^ = —V. This means that, although the gel thickness does 
not change, both interfaces are moving with the same velocity 
—V and therefore, v has the physical meaning of the tread- 
milling speed. 

If we now transform the dynamical equations into the co- 
moving frame moving with velocity in the direction of 
we can study the linear stability of the network thickness a 
with respect to perturbations of the type cos (mA^) at the in- 
ternal and external interface, e''(A^) and e^(A^), respectively, 

Fig.[T4]shows the dispersion relation for the largest growth 
rate depending on the wavenumber m. One of two eigenval- 
ues is always negative (i.e. stable), whereas the other one can 
be positive (unstable) depending on surface tension. We did 
not find a threshold value for the filament layer number, be- 
yond which the gel becomes stable towards small perturba- 
tions independent of the surface tension. However, a higher 
surface tension can suppress instabilities for thin gels. Re- 
cent experiments have shown the occuiTence of higher modes 
than one, i.e. the formation of up to three actin comets around 
one bead ifTj]], depending on the experimental conditions. 
The actual value for the surface tension of an actin network 
against water should be rather small, since actin is a solu- 
ble protein. Furthermore, typically small beads with higher 



IV. CONCLUSIONS AND PERSPECTIVES 

In this chapter we have tried to summarize the complex 
properties and out-of-equilibrium phenomena of actin gels 
linked by the Arp2/3 complex, which are at the origin of 
the motility of animal cells, as well as of intracellular or- 
ganelles and pathogens. Primarily, we have focused on two 
subjects; the complex actin polymerization dynamics under 
load at the polymer brush, and the symmetry -breaking of actin 
gels grown from the surface of small objects. While we have 
treated both subjects separately, it is obvious that a full under- 
standing of the system will have to include both approaches: 
the macroscopic stress distribution in the actin gel couples to 
the polymerization kinetics in the polymer brush, which in 
turn changes the deformation state of the gel and the macro- 
scopic stress distribution. We have shown in the previous 
paragraphs that homogenization models are at the moment 
the most appropriate models to capture the complex micro- 
scopic structures of biological materials on the one hand and 
take advantage of a continuous framework on the other hand. 
We believe, that the future in the modeling of growing actin 
gels in complex geometries, e.g. the advancing cell edge, lies 
in the coupling of these homogenization models to a com- 
plex dynamics in the polymer brush, as proposed e.g. in Refs. 

Ill ill]. 

An important question that remains to be elucidated is, how 
motion can be generated once the gel layer has become asym- 
metric. So far we have either limited our considerations to the 
case of a symmetry-breaking around objects, and neglected 
the generation of motion, or, as in most microscopic mod- 
els, we have only considered a stationary actin comet, which 
pushes an obstacle by polymerization. Both concepts are cir- 
cumventing, by more or less hand waving arguments, one crit- 
ical question. What is the origin of motion in the absence of 
external forces, provided that the actin comet and the object 
are only surrounded by a newtonian viscous fluid (recall that 
we are in a regime with Re ^ 1) and not attached to some 
support. Recently, Prost et al. ll48ll has put forward a simple 
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argument based on largely disparate friction coefficients for 
the obstacle and the actin tail and the property of treadmilling. 
In the following we will outline this argument. 

Suppose that an obstacle and its associated actin comet 
move with velocities Vo and Vc, respectively, in the labora- 
tory frame. In the viscous regime the force balance reads then 

= CoVo + CcVc , (45) 

where and denote the friction coefficients of the object 
and the comet with the surrounding fluid, respectively. The 
difference in the two velocities Vj = Vq — is the tread- 
milling speed. Substituting Vc in Eq. i45[ and after some re- 
arrangement one finds for the object velocity 

vo-vt-^. (46) 

t,o + Cc 

This means, that although the dissipation due to fluid friction 
is very small, it plays nevertheless a decisive role, since it is 
the ratio of friction coefficients which determines the object 
velocity. In the limit of ^ Co this velocity approaches 
the treadmilling speed. Given the fact, that the actin comet 
is much larger than the object, e.g. bead, droplet or vesicle. 



which causes a much larger friction, the experimentally ob- 
served obstacle velocities are indeed close to the treadmilling 
speed. 

Another biological aspect, which might be of importance 
when considering more complex cellular systems is the fact 
that the actin polymerization system is not constitutively 
active as in in vitro assays but is regulated by signaling 
cascades, which constitute in itself a nonlinear dynamical 
system. Typically these signaling cascades are modeled by 
reaction-diffusion systems which lead to pattern formation 
IstIi . e.g. polarization of the cell into leading and advancing 
edge. It remains to be shown how these two important 
mechanisms, elastic instabilities and instabilities due to 
reaction-diffusion processes, integrate to produce cellular 
motion. 
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